Crystal Structures of Polymerized Fullerides ACeo, A=K, Rb, Cs and Alkali-mediated 
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Starting from a model of rigid interacting Ceo polymer chains on an orthorhombic lattice, we 
study the mutual orientation of the chains and the stability of the crystalline structures Pmnn 
and 12 I'm. We take into account i) van der Waals interactions and electric quadrupole interactions 
between Ceo monomers on different chains as well as ii) interactions of the monomers with the 
surrounding alkali atoms. The direct interactions i) always lead to an antiferrorotational structure 
Pmnn with alternate orientation of the Ceo chains in planes (001). The interactions ii) with the 
alkalis consist of two parts: translation-rotation (TR) coupling where the orientations of the chains 
interact with displacements of the alkalis, and quadrupolar electronic polarizability (ep) coupling, 
where the electric quadrupoles on the Ceo monomers interact with induced quadrupoles due to 
excited electronic d states of the alkalis. Both interactions ii) lead to an effective orientation- 
orientation interaction between the Ceo chains and always favor the ferrorotational structure 12 jm 
where Ceo chains have a same orientation. The structures Pmnn for KCeo and 12 jm for Rb- and 
CsCeo are the result of a competition between the direct interaction i) and the alkali-mediated 
interactions ii). In Rb- and CsCeo the latter are found to be dominant, the preponderant role being 
played by the quadrupolar electronic polarizability of the alkali ions. 



I. INTRODUCTION 

Alkali metal doped Ceo (AajCeo), A=K, Rb, Cs, forms 
stable crystalline phases (fullerides) with a broad range 
of physical and chemical properties comprising supercon- 
ductors and polymer phases. For a review, see Refs. 
|^|-^. In particular the x = \ compounds Q exhibit 
plastic crystalline phases with cubic rock-salt structure 
(space group FmSm) at high temperature (T > 350 K) 
and polymer phases [||-^ at lower T. In the latter case 
the Ceo molecules are linked through a [2+2] cycload- 
dition 0, a mechanism originally proposed for photoin- 
duced polymerization |9| in pristine Ceo. From X-ray 
powder diffraction it was concluded that the crys- 
tal structure of both KCeo a-nd RbCeo was orthorhombic 
(space group Pmnn). Polymerization occurs along the 
orthorhombic a axis (the former cubic [110] direction). 
The orientation of the polymer chain, taken as a rigid 
unit, is characterized by the angle ■0 of the planes of cy- 
cloaddition with the orthorhombic c axis. In the Pmnn 
structure (Fig. la) the chains have alternating orienta- 
tions +1^ and —-0, 101 ~ 45 °. Notwithstanding this 
apparent structural similarity the electronic and mag- 
netic properties of KCgo on one hand and RbCeo and 
CsCeo on the other hand were found to be very different 
ESR and optical conductivity data [||,|o) show that 
RbCeo and CsCeo exhibit a transition from a quasi-one- 
dimensional metal to an insulating magnetic state near 
50 K, while KCgo stays metallic and nonmagnetic at low 
T. NMR spectra also did show marked differences be- 
tween K- and Rb-, CsCgo polymers [^. The contradic- 
tion between similar crystalline structures and different 
electromagnetic properties was resolved by single crys- 
tal X-ray diffraction and diffuse scattering studies [|2|. 



While the space group Pmnn is confirmed for KCeo, it is 
found that RbCeo is body-centered monoclinic with space 
group 12 /m. In the latter structure, the polymer chains 
have the same orientation 0^ (Fig. lb) in successive (001) 
pla nes. High-resolution synchrotron powder diffraction 
|l3| results have demonstrated that CsCeo has the same 
structure as RbCeo- Most recently a detailed structure of 
the polymer phase of K- and RbCeo has been performed 
by high-resolution neutron powder diffraction [ p^ . The 
distinct crystalline structures are confirmed and in ad- 
dition a determination of the positions of the C nuclei 
demonstrates a distortion of the Ceo monomers. The dis- 
covery of a metal-insulator phase transition by ESR spec- 
troscopy in KCeo around T = 50 K and the concomitant 
appearance of a superstructure revealed by X-ray diffrac- 
tion [|l^] demonstrate a subtle interplay of structure, di- 
mensionality and electronic properties. In fact the nature 
of electrical conductivity in the ACeo compounds is still 
under debate. Previous electronic band structure calcu- 
lations (lelji^l suggest a three-dimensional dispersion of 
the electronic bands. However numerous experiments 
are interpreted in terms of a quasi-one-dimensional con- 
ductor. 

Concerning our understanding of the crystalline struc- 
ture, a condensation scheme has been proposed for 
the phase transition from the orientationally disordered 
structure Fm3m to the orientationally ordered polymer 
phase Pmnn on basis of a phenomenological Landau 
theory [|l^. However, a microscopic description of this 
transition (which involves the cycloaddition process) is 
very complex. Already in the FmSm phase, the charge 
transfer of one electron from the alkali atoms to the Ceo 
molecules giving rise to a A+C^g crystal changes the elec- 
tronic structure of the molecule. This change affects the 
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FIG. 1. Crystal structures projected onto the crystallo- 
graphic {b,c) plane: (a) Pmnn, (b) 72/m. The thick bars 
represent the projection of the cycloaddition planes. Poly- 
merization occurs along a. The alkalis located in (6, c) planes 
and at ±a/2 are denoted by full (+) and empty (— ) circles. 



orientation dependent interactions, in particular the crys- 
tal field, and favors polymerization The latter is 
a quantum-chemical process and within our view it is 
the driving mechanism for the phase transition. There 
remains the question why polymerized KCgo and Rb-, 
CsCgo have different crystal structures and different elec- 
tronic properties and why pressure polymerized Cgo has 
the same crystalline structure as KCgo. In the follow- 
ing we will essentially address the problem of the various 
crystalline structures. It turns out that the differences 
in structure can be related to the average radii of ex- 
cited electronic c?-states of the alkali cations. A prelim- 
inary account has been given as a letter |2^. There it 
was found that the direct electric quadrupole interactions 
between charged polymer chains in the alkali fuUerides 
(the existence of electric quadrupoles on the chains is 
a consequence of the charge transfer between an alkali 
atom and a Cgo monomer [ ^ ) always favors the struc- 
ture Pmnn with alternating chain orientations. While 
this result agrees with experiment for the case of KCeo, it 
does not explain the monoclinic structure 12 jm of RbCgo 
and CsCgQ. Hence it becomes necessary to take into ac- 
count the role of the alkalis. In particular the specific 
quadrupolar polarizability of the alkali metal ions leads 
to an indirect interchain coupling which favors 11 jm with 
equal chain orientation. The competition between direct 
and indirect interactions then explains the structural dif- 
ferences between KCgo and RbCgo, CsCeo- In the present 
paper we will present several important extensions of the 
previous work. 

Since polymerization leads to chains of D2h symme- 
try, one has to investigate the role of van der Waals in- 
teractions (repulsive Born-Mayer and attractive London 
dispersion forces) for a configuration that is very differ- 
ent from the situation in orientationally disordered Cgo- 
Potential energy calculations based on van der Waals 
forces between Cgo polymer chains show that the energy 



minimum is highly sensitive to the relative chain orien- 
tation We will show that the dominant term of 
the multipole expansion of the van der Waals interaction 
between different chains, corresponding to quadrupole- 
quadrupole interactions, always favors the orthorhombic 
Pmnn structure too. Our study of the van der Waals 
interchain interactions is not only relevant for the struc- 
ture of the ACgo compounds, where the direct interchain 
interactions are in competition with the alkali-mediated 
chain-chain interactions, but also for the understanding 
of the structure of polymerized Cgo . There a "low pres- 
sure" orthorhombic phase was found to have the struc- 
ture Pmnn [|3|j24). 

As a further extension we will study the coupling be- 
tween rotational motion of the polymer chains and the 
lattice translations of the alkalis (TR coupling) . Thereby 
we will obtain the shear mode that characterizes the mon- 
oclinic unit cell of RbCgo @ and CsCgo (iS). How- 
ever our calculations show that the indirect chain-chain 
interaction mediated by the displacements of the alka- 
lis is relatively weak and is not able to account for the 
structural difference between KCeo and Rb-, CsCgo. The 
quadrupolar polarizability mediated interaction is the de- 
cisive mechanism. 

The content of the paper is as follows. In Sect. [| we 
formulate the rigid chain model of Cgo polymer chains 
in the orthorhombic lattice (space group Immm). The 
sole degree of freedom of each chain is the rotation an- 
gle ip about the long axis. Since the interchain poten- 
tial is a function of the angles tp, we perform a multi- 
pole expansion into symmetry adapted rotator functions 
(SARFs), taking into account the symmetry of the poly- 
mer chain (-D2/1) and of the site. The direct interchain 
potential is studied for both Coulomb and van der Waals 
forces between monomers on different chains (Sects. || 
and IH). These forces lead to Coulomb- and van der 
Waals type quadrupole-quadrupole interactions between 
chains. Studying the direct interaction potentials in 
Fourier space, we find that they always favor the struc- 
ture Pmnn. Next (Sect. IV) we study the interactions of 
the Cgo chains with the surrounding alkalis. We consider 
two types of interactions. Firstly the coupling of the rota- 
tions of the polymers with the alkali displacements (TR- 
coupling) via van der Waals and Coulomb forces, and 
secondly the coupling of rotations of the polymers with 
induced quadrupoles on the alkali cations (quadrupolar 
polarizability coupling). Both interactions lead to an ef- 
fective rotation coupling between polymer chains, which 
always favors a ferrorotational structure /2/m. In the 
following, we study the competition between direct and 
indirect interactions and discuss the stability of struc- 
tures (Sect. Finally (VI) conclusions are drawn. We 
discuss the unique properties of ACgo in comparison with 
ionic molecular crystals with small ions. The paper has 
three appendices, where we discuss: A) the orthorhom- 
bic lattice structure as a consequence of polymerization, 
B) details about the interchain potentials, C) the mi- 
croscopic origin of the quadrupolar polarizability of the 
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FIG. 2. Fragment of a Ceo polymer chain, orientation 

= 0. 



cations. 



II. RIGID CHAIN MODEL 



We start from the assumption that the polymeriza- 
tion by stereospecific cycloaddition of Cgo molecules in 
ACeo has occurred along the original cubic [110] direc- 
tion. Polymerization acts as a negative uniaxial stress 
along [110]. By using concepts of the theory of elasticity 
psi , we find that the cubic crystal is deformed into an 
orthorhombic one (point group D2h)- Taking the cubic 
[110], [110], and [001] as new ^, ry, C, axes (orthorhombic 
a, h, c) respectively, we find the deformations 
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[cii (cii -I- ci2 -I- 2cm) - 2C12] , 

[cil (Cii -I- C12 - 2C44) 



(2.1a) 

(2.1b) 
(2.1c) 



where d = cn (cn -I- C12) — 2cf2, K > 0. Details of the 
derivation are given in Appendix From the relations 
between the cubic elastic constants we see that e^j < 0, 
Criri > and e^^ > 0, which corresponds to a contrac- 
tion along the ^ direction (which we identify with the 
orthorhombic a axis) and to elongations along the rj and 
C directions (the orthorhombic b and c axes respectively). 
In the following we consider an orthorhombic lattice with 
space group Immm, with polymer chains oriented along 
the a axis (Fig. 2). We take the chains as rigid units of 
D2h symmetry where the sole degree of freedom is the ro- 
tation angle "0 about the chain axis a. The assumption of 
a rigid chain is a reasonable first approximation, indeed 
vibrational density of states data on RbCgo obtained by 
inelastic neutron scattering exhibit a low energy external 
mode below 5 meV which is interpreted as arising from 
librations around the chain axis [ p6[ . Inelastic neutron 
scattering results on the orthorhombic phase of KCgo are 
in close agreement with these results [ pTj ]. 

The mathematical formulation of our model of a poly- 
merized crystal is as follows. We consider N Cgo 
molecules located on a body-centered orthorhombic lat- 
tice, where the center of mass positions of the molecules 



coincide with the lattice points. The lattice points are la- 
beled by indices ft = (ni, 712, na). These are either all in- 
tegers or all half integer numbers, corresponding respec- 
tively to the corners and the center of the orthorhombic 
cell. The equilibrium position of the n-th lattice point 
then reads 



X{n) — nid + 7126 -|- n^c, 



(2.2) 



where a, b and c are the orthorhombic lattice vectors. 
The molecules (monomers) in a polymer chain are ori- 
ented so that a twofold molecular axis coincides with the 
orthorhombic a axis (^-direction). The number of chains 
is Nc, and the number of monomers within a chain is iVi, 
hence N = NcNi . We define as standard orientation of a 
polymer chain the orientation where the plane of cycload- 
dition is parallel to the (a, c) planes of the orthorhombic 
crystal. The angle ip then measures a counterclockwise 
rotation of the polymer about a, the standard orienta- 
tion corresponding to V' = (Fig. 2). Since it is sufficient 
to specify the indices v = (n2,n3) to address a polymer 
chain, and because all molecules along one chain have the 
same orientation, the rotation angle is independent of the 
index ni and is written as ip = i/'(j7). The interaction be- 
tween chains is then formulated as a sum of two-body 
potentials U{n,ip(iy);n',ip{i'')) between Cgo monomers: 



(2.3) 



This interaction depends on the mutual orientation of 
the chains. In order to describe the symmetry reduction 
(phase transition) from a hypothetical polymer crystal 
with structure Immm (i. e. no preferential orientation 
angle ipii^)) to Pmnn or /2/m, we perform a multipole 
expansion of the potential in terms of SARFs Si{i') = 
5z(V'(z7)) = sin(;V(z7)), Q(i7) = Q(^(z7)) = cos{l^{,7)), I 
being the angular momentum quantum number. SARFs 
are the most appropriate variables for the description 
of orientation dependent properties in molecular solids 
[^1^. Since in the present problem the assumption of 
a rigid chain leads to one single rotation angle our 
SARFs are particularly simple and correspond to pla- 
nar rotator functions [ |3l| ]. Symmetry of the chain for 
■0 — > ip + t: implies that only even values of I occur. In 
lowest order of the multipolar index I we obtain 



RB. 



VcF, 



(2.4) 



with 



VqQ = ^ E E JQQi^^^^-^ i7')^2(j7)52(?'), (2.5a) 

^cF^Y^Y. <^^')C2{iy). (2.5b) 

Here Vqq stands for the quadrupolar (/ — 2) rotation- 
rotation interaction The quadrupolar function 5*2, which 
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is uneven in ipi^v), is an appropriate order parameter. 
The coefficients Jqq are obtained from 

jQQ{ni,ty;n[,i/) ^ — dip{v) dip{v') 
I" Ja Ja 

X U{n,ij{i7);n\ib{P'))S2{i^{iy))S2WP')). (2.6) 

The term Vcf accounts for the crystal field. It is the 
potential experienced by a central chain with orientation 
-(/^(i?) where the surrounding chains are taken in cylindri- 
cal approximation {I = 0). The coefficients v are given 
by 

X C/(n,V(z7);n',V(^7'))C2(^(^7)). (2.7) 

For a Ceo monomer on a site (ni, u) we consider the in- 
teractions J a with the two monomers at {jii ± i) on the 
four neighboring chains (n2 ± ^,n3 ± ^) and the inter- 
actions Jf, with the one monomer on each of the chains 
(n2 ± 1,713). 
We introduce Fourier transforms 
1 



S2{q) 



,iq-X{n) 



(2.8a) 
(2.8b) 



Since the chains are rigid and form a two-dimensional 
array v, S2{n) = S2{ni,i') is independent of ni, and 
hence S2{q) reduces to 



— ^52(n)e-^-^("). 



52(0, gl) = v^<5,,o52(gl) 



with 



S2{qi 



E^2(?) 



(2.9) 
(2.10) 



Here we have X{iy) — + n^c, and q± = (^rp^c)' 94' 
q^, (j^ being the components of q along the orthorhombic 
axes. The interaction then reads 



RR 
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E 



jQQ{q±)Sl{ql)S2{qi 



^QQ('7-L) = 8JaCos('^jcos(^) +2JbCos(g„6). (2.12; 



(2.11; 



In the following section we will derive explicit expres- 
sions for the quadrupolar interaction coefficients Ja, Jt 
and for the crystal field. The crystal field for the same 
configuration of neighboring interactions reads 



VcF^Niy(8va + 2vi,)C2ii^), 



(2.13) 



where the coefficients Va and Vh are determined by Eq. 
(2.7). In the next section the theory will be applied to 
the study of van der Waals and electrostatic quadrupole 
interactions between Cgo polymer chains. 



III. DIRECT INTERCHAIN COUPLING 

The orientation dependent interaction between two 
polymer chains in the Immm crystal of ACeo is due to di- 
rect van der Waals interactions and to electric quadrupole 
forces between Ceo monomers. The van der Waals inter- 
actions also determine the mutual chain orientations in 
polymerized Ceo- 



A. Direct van der Waals Interactions 

Here we consider repulsive Born-Mayer and attractive 
London dispersion forces between interaction centers on 
monomers located in different chains. The interaction 
centers comprise C atoms and double bond- and sin- 
gle bond centers; they have been previously introduced 
for the description of the orientational disordered and 
ordered phases of Ceo fuUerite While in ori- 

entationally disordered Ceo sll C atoms are equivalent, 
this is no longer the case in the polymer phase and we 
have to take into account explicitly the geometrical con- 
straints imposed by polymerization along the orthorhom- 
bic a axis. Throughout this paper we assume polymer 
chains of symmetry D2h- The neutron diffraction results 
of Ref. |l4j show in fact that the symmetry of the chains 
in both KCeo and RbCgo is lower, C2h- We attribute this 
reduction of symmetry to the molecular field of the origi- 
nally ordered structures Pmnn and 12 /m. Since we start 
from an orientationally disordered structure Immm and 
study the transition towards Pmnn or 72/to, we keep 
the chain symmetry D2h- It is convenient to classify the 
C atoms of a monomer in planar (100) sheets, perpendic- 
ular to a (see Fig. 1). We label these sheets by an index 
A = 1, . . . , 17. The position of sheet A with respect to the 
center of mass of the monomer is determined by a vector 

pixy. 



p{X) =p{X)e^, 



(3.1) 



with origin at the center of the monomer and end point 
at the intersection of sheet A with the polymer axis. This 
intersection is called the center of the sheet. Here is 
a unit vector along the polymer axis. The atoms within 
sheet A are labeled by an index c(A). We will use the 
notation A = (A,c(A)). The atomic positions relative to 
the center of mass of a monomer belonging to the polymer 
iJ which is rotated away form the standard orientation by 
the angle then reads 



f(A, ^P{P)) ^ r{X){cos[(t>{A) + V^(j7)]e^ 
+ sin[(/)(A)+i/;(i7)]ec} 



(3.2) 



with r(A) = d sin 0(A). Here d is the radius of the Ceo 
molecule, 6{X) and (j){A) are respectively the polar and 
the azimuthal angle of the atom c(A) when the polymer 
is in the standard orientation. The position of atom A of 
the monomer at lattice site n — (ni,i') is then given by 
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TABLE I. Lattice constants of the cubic and orthorhombic 
lattices of ACeo and Ceo (units A). Calculated direct interac- 
tion coefficients J for van der Waals and Coulomb interactions 
(units Kelvin), orthorhombic lattices. 



TABLE IL Crystal field coefficients due to neighboring 
chains v}^ etc., units K. 





flc ao bo 


Co 


tW" tW 
T VV T VV 

■Ja ''h 


Ja 




KCeo 


14.06 9.11 9.95 


14.32 


21.44 -6.61 


10.70 


-55.17 


RbCso 


14.08 9.14 10.11 


14.23 


21.33 -2.88 


11.14 


-50.56 


CsCeo 


14.13 9.10 10.22 


14.17 


20.61 -1.50 


11.65 


-47.68 


Ceo 


14.15 9.14 9.90 


14.66 


13.13 -8.41 


/ 


/ 




i?(n,A,V(i7)) = 


Xin) 


+ p(A,?/.(j7)), 




(3.3) 
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w 
















KCeo 


24.24 


-73.06 


-99.12 


422.53 


RbCeo 


23.44 


-15.31 


-92.99 


400.99 


CsCso 


21.49 


5.57 


-89.55 


387.06 


Ceo 


24.45 


-101.01 


/ 


/ 



where p{A,ip{i')) — p(A) + r{A,ip{i')). The distance be- 
tween atoms A and A' belonging to monomers ft and n' 
in different chains i7 and reads 

A(n,A,V;n',A',V') = \R{n,A,tp)~ R{n',A',il;')\. (3.4) 

Here we have written ■0, V'' for "fpi^) ^^nd ipi^') respec- 
tively. The van der Waals potential between these 
atoms is then given by 

U^{n, A, ij; n', A', V') = Ci exp(-C2A) - BA'^ (3.5) 



with A given by Eq. (3.4). Here Ci and C2 are the pa- 
rameters of the repulsive Born-Mayer potential while B 
determines the strength of the van der Waals attraction 
(London dispersion forces). The interaction potential be- 
tween monomers then reads 



B. Electrostatic Quadrupolar Interactions 

The charge transfer of one electron form the alkali 
atom to the Cgo molecule leads to an occupation of 
the lowest unoccupied molecular orbital (LUMO) levels 
which are of Tin symmetry ]3^ . Thereby the crystal field 
of the C^o ion acquires an electronic component that fa- 
vors the same orientation of the neighboring molecules 
along [110] such that the stereospecific cycloaddition oc- 
curs [H. We have studied the electronic charge distri- 
bution on the Cqq units in the polymer chain by using a 
tight binding model pl[ |. The charge is mainly concen- 
trated in the equatorial region of CgQ, in agreement with 
NMR results |3^. We find that only even / multipoles 
are allowed; in particular each Cgo unit has an electric 
quadrupole. In the following we adopt a simple model of 
charge distribution. By using the labeling of C atoms of 
we locate a charge of —0.12 (units electron charge 
|e|) on each bond C15-C16. These charges are fixed at a 
distance d = 3.52 A from the center of the C^g ball. In 
the center we put a charge —0.76. The position of these 



U^in,^;n',i;') = Y,U^'iri,A,'ilj;ri',A','tP'), (3.6) three char ges on a CgQ monomer centered at the lattice 



A, A' 



where the double sum runs over the C atoms A, A' of 
the two monomers. In a similar way we can treat the 
potential due to interaction centers located on double 
bonds and single bonds of the monomers. For the case 
of three interaction centers between two C atoms partic- 
ipating in a double bond on the monomer, the present 
formulation leads to the consideration of 26 additional 
planar sheets, the inclusion of single bond centers adds 
17 more sheets. Details of the van der Waals potentials 
between various types of interaction centers belonging 
to monomers on different chains are given in Appendix 
As a result the potential between two monomers is 
again expressed by Eq. ( |3.6[ ), where now the indices A 
and A' run over atoms, double bond- and single bond 
interaction centers. Having dete rmin e d U^ (n, ip; n' , ip'), 
we proceed as in Sect, g Eqs. (p3|)-(p^). As a result 



site fi — [rii , v) belonging to chain v is then given by 



i?(n, a, V(!7)) = X{n) + D{a, 7/;(z7)), 



(3.7) 



where 



D{a, tpiv)) = d{a) [e^ cos^]{v) + sin ^'(i^)] , (3.8) 

where d{a) = ±d for a — 1,2 corresponding to the two 
charges Qa = —0.12 and d{a) = for a — 3, correspond- 
ing to the central charge Qa = —0.76. The distance 
between two charges belonging to two monomers on dif- 
ferent chains is then given by 

A(n, a, ■0; n , a\ ifj') — \R{n, a, -0) — i?(n', a', V'OL (3-9) 

where V'j "0' stands for V'C^); The Coulomb inter- 

action between these two monomers then reads 



we obtain the van der Waals contributions and 
to the quadrupolar interactions J a and Jfj. Numerical re- 
sults are found in Table |. There we quote the interchain 
van der Waals interaction coefficients for orthorhombic 
lattice constants Co, 60, Co of polymerized ACeo, A=K, 
Rb, Cs and polymerized Cgo- Similarly we determine 
the crystal field coefficients and ■ The results are 
quoted in Table ||. 



f/^(n, .A; n', 0') = V p-—^^ -. 

^ A(n, a, ■0; n', a', -0') 



(3.10) 



With our units of charges and with lengths in A, with 
F = 167000 K A, the energy U'~^ is measured in units 
Kelvin. One now proceeds as before with . We take 
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into account the same configuration of neighbors as be- 
fore, sum over all chains in the crystal, expand the to- 
tal potential in terms of SARFs and obtain the interac- 
tion coefficients , jf for the electrostatic quadrupole- 
quadrupole interaction and the coefficients and 
for the electrostatic crystal field. The results are also 
quoted in Table | and Table || respectively. Adding 
the quadrupole-quadrupole contributions from the van 
der Waals and the electrostatic interchain potentials we 
obtain the total quadrupolar interaction Vqq, expres- 



sion (2.11), where JQQ{qi_) is given by Eq. (p.l2|) with 



J a = 
Jb 



jW 
-'a " 



J? 



a ? 
■C 
b ■ 



(3.11a) 
(3.11b) 



Notice that for polymerized pristine Cgo, there is 
no electrostatic quadrupole due to charge transfer, and 
hence and are zero. We observe that in principle 
the deformation of the monomer in Ceo from — » D2h 
leads to a redistribution of the electric charge on the neu- 
tral molecule. Here it is assumed that these effects are 
included in and . 

We return to the quadrupolar interaction J{q) and 
study its wave number dependence. In two-dimensional 
q± space (g^, g^) we have at the Brillouin zone (BZ) cen- 
ter q± = (0, 0) = gr 



jQQiqr)=SJa + 2Jb 



(3.12) 



With the numerical values of Table | we conclude that for 
polymerized ACgo, A=K, Rb, Cs as well as for polymer- 
ized Ceo, Jiqr) > 0- On the other hand, at the Brillouin 
zone boundary q± ~ {0,^) — qz we get 



Jqqiqz) — —8Ja + 2Jf,, 



(3.13) 



with J{qz) < and |J((f2)| > J{qr)- The dominance 
and negative sign of J{qz) leads to a condensation of 
S2{q) at ql = qz'. 



{S2{q±)) = VNcaSif^^if^ 



(3.14) 



which is the condensation scheme proposed earlier [[l8[ . 
Here ( ) stands for a thermal average over the crystal, a 
is the order parameter amplitude. The number of chains 
Nc corresponds to the number of lattice points in the 
orthorhombic (6, cj plane. Condensation at qz implies 
that the chains in the same basal plane (a, b) of the or- 
thorhombic lattice all have the same orientation, but the 
orientation alternates in neighboring planes at distance 
c/2. This is the "antiferrorotational" structure Pmnn 
(Fig. la). We find that the rotation-rotation interaction 
due to van der Waals forces and electrostatic quadrupolar 
forces leads to Pmnn for KCeo, RbCeo and CsCeo, irre- 
spective of the orthorhombic lattice constants. While for 
KCgo this structure has been found by experiment ||7|, |l2| , 
the experimental result for RbCgo jl^jl^l CsCgo |l3| ] 
is 12 /m. We conclude that the alkalis must play a spe- 
cific role and we will investigate this problem in the next 
section. 



On the other hand, for polymerized pristine Cgo, 
our analysis confirms the experimental result Pmnn 
11,11. 



IV. ALKALI-MEDIATED INTERACTIONS 

A. Translation-Rotation Coupling 

In orientationally disordered molecular crystals or in 
crystals containing ions with orientational disorder, the 
coupling of molecular rotations with center of mass lat- 
tice displacements of molecules and of surrounding ions 
plays an important role in determining the structural 
properties. For a review on the so-called translation- 
rotation (TR) coupling, see Ref. Here we will study 
the corresponding coupling of the rotational motion of 
polymer chains with the lattice displacements of the sur- 
rounding alkali ions. We extend the model of an or- 
thorhombic crystal consisting of orientationally disor- 
dered rigid polymers with axis of polymerization along 
the orthorhombic ^ direction (Sect. ^) by including the 
alkali metal ions at equilibrium lattice positions 



X{nA) = nAia + nA2b + nA3C, 



(4.1) 



where Ha = (j^Ai, ?^A2j ^^As) are the lattice indices of 
the alkali atoms. For a CgQ monomer at lattice site 
n = (ni,n2,n3), the six nearest neighbor alkalis are lo- 
cated at (ni ± i,n2 ± ^jUs), (ni,n2,n3 ± i). Taking 
into account lattice displacements ma(?ia) we write for 
the actual positions of the alkalis 



R{nA) = X{nA) + UAinA)- 



(4.2) 



The distance between a C atom in position A belonging to 
a CgQ monomer at lattice site n = (rii, j7) in the polymer 
chain i7 and a surrounding alkali ion at Ha is given by 
\R{n,A,ip{H)) — R{nA)\ or equivalently by 

A(n, A,'(/'(i');nA) = \p{A,i>{v)) + X{n~nA) - UA{nA)\. 

(4.3) 



Here we have used expression ( |3.3[ ) for R{n, K,jp{v)), 
TpiH) being the orientation angle of the polymer chain. 
The van der Waals interaction of the monomer and the 
alkali ion is then given by 



U'^in,iPiH);nA) 



A 



Ci exp(-C2A) - BA-^, (4.4) 



where the potential parameters Ci, C2 and B refer to 
the C-A interaction (see Appendix The Coulomb 
interaction potential of the quadrup olar e lectric charges 
Qa on the C^q monomers (See Sect. IIIB) with an alkali 
ion A+ with unit charge reads: 



U^in,i,iH);nA)^ fJ2 



QctQ 



A{n,a,A,ijj{iy);nA) 



(4.5a) 
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where 



= \D{i,ip{v)) + X{n- Ha) - ma ("a)! 



(4.5b) 



is the distance between the charge a on the monomer 
and the center of the alkah ion. Assuming that the lat- 
tice displacements are small (in comparison to the lattice 
constants), we expand these potentials {U stands for 
or U^) in terms of ua and retain the first two terms: 



U{n,ij{iy);nA) = U^"Hrt-rtA,^{iy)) 
+ ^U^^\n~ nA,'4'{v))uAt{nA), 



(4.6) 



where i labels the Cartesian components along the or- 
thorhombic axes. The first term on the right hand side 
(RHS) corresponds to the situation ua = 0, the alkali ion 
is at its equilibrium position. Exploiting the angular de- 
pendence of /7^°-' we expand it in terms of SARFs. Since 
each monomer is located at a center of symmetry with 
respect to the surrounding alkalis in the orthorhombic 
lattice, only Ci{;i]-){v)) terms contribute to 

Y,U^°\n-nA,^{v)) = ^z;(°)(nA - n)C2{^{v)), 

Ha nA 

(4.7) 



where 

v^^^ (riA — ft) 



d^(i;){/(0)(7?-n^,V;(i7))C2(^(i7)). 

(4.8) 



Here we restrict ourselv es t o the lowest order multipoles 
with 1 — 2. Expression (4/7) stands for the interaction of 
one monomer at site ni in chain ly with six surrounding 
alkali atoms at equilibrium lattice sites X(n^). Since all 
monomers of the chain are equivalent, the RHS result 
does not depend on ni. Taking into account U'~^ and 
and summing over all monomer sites, we obtain the 
crystal field potential due to the nearest neighbor alkali 
atoms 



Vc^ ^ N,J2i^v^t 



2vf)C,mV)) 



where we have 



A,W 



A,C 
ab ' 



(4.9) 



(4.10a) 
(4.10b) 



Here Vab and Vc refer to alkalis located at (±i, ±i, 0) and 
(0, 0, ±i) respectively, counted from the monomer center 
n. Numerical values of these coefficients are quoted in 
Table Ell 



TABLE III. Crystal field coefficients due to alkalis, units 



K. 





A,W 




A,C 




KCeo 


-10.55 


16.38 


-497.98 


1137.42 


RbCeo 


-14.30 


26.44 


-493.28 


1161.03 


CsCeo 


-19.91 


30.44 


-497.81 


1177.16 



The second term on the RHS of Eq. (4.6) is the proper 
TR contribution, u\^^ stands for the first order deriva- 
tive of U with respect to UAi- Expanding ul^^ in terms of 
SARFs, summing over the surrounding six alkalis of the 
monomer and exploiting the symmetry of the orthorhom- 
bic lattice, we find in lowest order ^ — 2 of the multipole 
expansion: 



'S^u'l^\n- nA,ii{v))uAi{nA) 



UA 



^w|^^(riA - n)S2{ilj{v))uAi{nA)- 
Ha 



(4.11) 



Here we retain only the terms containing the function 
^2 and drop those with C2. Indeed it can be shown 
that TR terms with C2 do not contribute to a change 
of the orthorhombic lattice structure. The coefficients 
i^i^^T^A ~ n) are obtained by 



.,(1) 



[nA - ft) = - 

TT 



d'^{i7)U^'\n~nA,ij{i^))S2{^j{i^)). 

(4.12) 



Summing over all monomers in the crystal, the total TR 
potential due to + U'^ is found to be 



V 



TR 



= EEE^» (^^ - n)S2im)nA.iriA). (4.13) 



n riA I 

We mention the symmetry property 

(1)/- -^ (1)^- - ^ 

vl '[UA -n) = -v] '(n - ua). 



(4.14) 



From expression ( 4.12| ) it follows that all six coefficients 
v'^^ (nA — n) are zero. The only nonzero coefficients 
are v^i\o,0,±^) = T<c, "^'^(±i±iO) 



„(i) 



(±i,TiO) = T 



"3,ab 



3,a6' 
A,W 



3 V^2' 

where 



V. 



V. 



2,c 
.A,W 



( A,C 



A,C 



Numerical values are given in Table IV 
Fourier transforms of displacements by 



UA{nA) 



UA{q) 



1 



XiuA) 



UA{q), 




UA{nA)- 



Tvi,ab and 



(4.15a) 
(4.15b) 
We define 

(4.16a) 

(4.16b) 
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TABLE IV. Calculated indirect interactions coefficients 
for TR coupling [v in units K/A). Electronic polarizability 
parameters, dA in units A, qa and calculated A;,, Ac in units 
K. 





A,W 


A,W 


A,C 


A,C 
''^S.ab 


d-A 


At 


Xc 


9A 


KCeo 


4.57 


4.24 


317.71 


199.98 


1.47 


9.28 


51.20 


43.25 


RbCeo 


7.48 


5.70 


326.36 


195.17 


1.82 


13.25 


83.52 


35.00 


CsCeo 


8.59 


7.79 


332.29 


194.84 


1.87 


14.27 


90.75 


34.00 



Using in addition Eq. ( 2.10 ), we rewrite Eq. ( 4.13 ) 



TR 



V 



where 



J2 E vl'\q±)Sl{q^)uM{0, gl), (4.17) 



rriA '-Z 

q± « 



y^^\'^±) = *E^»^'^('"^-4 - n)sin(q-- X{nA ~ n)), (4.18) 



or equivalently 







'4w^abSin(9r,|) 



(4.19) 



We see that the rotational motion of the chains about 
the a axis induces displacements of the alkali atoms only 
along the b and c axes. We observe that the coupling 
vl^\q±) vanishes at q± — qz- 

In order to relate our results to elastic deformations of 



the lattice, we consider the RHS of Eq. (4.17) in the long 
wavelength regime and transform to acoustic (ac) lattice 
displacements 



s{q) = \ —UA{q) 



m 



(4.20) 



where m = niA + 'm-Ceo total mass per primitive 

unit cell, mcgo being the mass of a Ceo monomer. The 
acoustic part of V^^ then reads 

= VA^EE^^*'^(9^)^^(0''?'^)^2(9l). (4.21) 

The sum over i is now restricted to the components 77 
and C and hence 



(4.22) 



The translational part of the acoustic lattice energy is 
given by 



(4.23) 



q «J 



Here Mij{q) is the dynamical matrix of the orthorhom- 
bic crystal in absence of TR coupling. In the following 



TABLE V. Interactions at 



gr, units K. 









— Jac{qr) 


-Jep{qr) 


■Hqr) 


KCeo 


-24.80 


158.30 


-37.72 


-225.02 


-129.24 


RbCeo 


-12.03 


164.88 


-39.44 


-691.62 


-578.21 


CsCeo 


-2.18 


161.88 


-40.65 


-837.14 


-718.09 


Ceo 


/ 


88.22 


/ 


/ 


88.22 



we need only to retain the components (77, ^) of s(0, gl). 
Then M{q) reduces to a 2 x 2 matrix: 



M{qi_) = 



P 



qf,C22 



qlcii qr,qc{c23 

g»)<7c(c23 + C44) 



+ Ca4) 

q^cAA 



(4.24) 



Here p is the mass density in the body centered or- 
thorhombic unit cell, p = 2m/{ahc). We use the Voigt- 
notation for the orthorhombic elastic constants C22 etc.. 
We now consider 



Vnr = V™ + V 

* 0.C * ac ' ^ ac 



TT 



(4.25) 

For a fixed configuration of orientations {S2{qi_)}^ we 
minimize Vac with respect to Si(0,q±), i — rj^Q, and ob- 
tain 

st(0,9l) = -VA^Af-i(gl)«(i)(gl)52^(gl). (4.26) 

Substitution of this expression into Vac leads to the ef- 
fective rotation-rotation interaction which we denote by 

yRR. 



V 



RR 
ac 



9J- 



)S2{qA 



where 



Jac{qi_) = ^j^(?-l)Af-Hgl)eH'a). 



(4.27) 



(4.28) 



Since Jac{q±) > 0, the lattice mediated interaction V^^ 
is always attractive. The largest value is obtained for 

q±= qr = 0: 



^^Ihn^ Jac(gl)|,^=o 



caab 



Jaciqv)- (4.29) 



With 



C44 



870 K A- 



^ taken from Ref. 



for C, 



60" 



fuUerite, and the values of c from Table [V, we find 
Jaciqv) as quoted in Table ^ for KCeo, RbCeo and CsCeo. 
The attractive interaction at qr favors a condensation of 
a ferrorotational structure where all polymer chains have 
the same orientation, i. e. the space group /2/m. How- 
ever comparison of the numerical values of Jac{qv) and 
J{qz) shows that the interaction, which leads to the anti- 
ferrorotational structure Pmnn, is dominant. Hence the 
TR coupling mechanism is not sufficient to explain the 
stru ctural difference between KCeo and Cs-, RbCeo. In 
part [V B we will exploit a different alkali-mediated inter- 
action mechanism. It is based on a specific quadrupolar 
electronic polarizability of the alkali ions. 
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B. Quadrupolar Polarizability 

We will now include the role of the specific electronic 
polarizability of the alkali ions. It is known from work 
on the ammonium halides p9[ that the indirect interac- 
tion of two NH^ tetrahedra via the polarizable halide 
ions plays an essential role in determining the various 
crystalline phases of the ammonium halides NH4X, X = 
CI, Br and I. However, in the present problem, since the 
Cgo monomers have symmetry D2h, they do not couple 
to the dipolar electronic polarizability of the alkali metal 
ions. We have to resort to the quadrupolar polarizabil- 
ity. Since the C^q units in a polymer chain are rigidly 
linked in the same orientation, the CgQ chains, which sup- 
port electric quadrupoles, produce coherent electric field 
gradients which induce an anisotropic (quadrupolar) de- 
formation of the residual electronic charge of the alkalis. 
The presence of residual electronic charge is seen as a 
consequence of the uniquely large interstitial space be- 
tween the Ceo molecules (we recall that in ACgo, the A+ 
cations occupy the formerly octahedral positions of the 
cubic phase). Within a touching sphere model, with 5 A 
as the effective van der Waals radius of the Cgo molecule 
po| and a cubic lattice constant a = 14.15 A, we estimate 
the radius of the interstitial sphere to be Ra ~ 2.075 A. 
Although the charge transfer from the alkali atom to the 
Ceo molecule is considered as complete, there still will be 
a charge of 0.1 — 0.15|e| left inside the interstitial sphere 
centered at each alkali site. In Appendix |^ we discuss 
in more detail the microscopic origin of the quadrupo- 
lar polarizability of the cations. Within a tight-binding 
model of the conduction electron band, we expand the 
electron wave function at an alkali site in terms of lo- 
cal s and d functions. If there is an appreciable weight 
of d-states, then the alkalis acquire a quadrupolar mo- 
ment. We model the corresponding charge distribution 
of each alkali ion by a symmetric linear dumbbell cen- 
tered on lines along the a direction. We take dumbbells 
with equal charges = qa, /3 = 1, 2, at distances ±dA 
from the center. The numerical values of dA (Table IV) 



are the average radii of valence electron d shells calcu- 
lated with atomic wave functions 8^3/2, 4^3/2 and 6^3/2 
for K+, Rb"*" and Cs"*" respectively (Appendix ^). We 
observe that the dA values for Cs"*" and Rb+ are close to 
each other but differ from K+. We consider d shells be- 
cause they can support an electric quadrupole moment. 
On a same crystalline line along a, these dumbbells are 
parallel with their axis perpendicular to a and a same ori- 
entation angle tp with the c axis. In the crystal we then 
have chains of alkali dumbbells, where the rigid chain 
is not imposed by intrachain interactions (as is the case 
for the Cgo polymers formed by cycloaddition) but by 
the surrounding C^q chains. The location of a dumbbell 
in the crystal is determined by Ha — ('t-aii^'a), where 
ua = (?^2A,"'3a) denotes the chain and where uai la- 
bels the dumbbell within the chain. The orientational 
motion of an alkali dumbbell is characterized by SARFs 



•S2(^'a) = sm{2'ip{i7A)), independent of riAi- The electric 
quadrupole-quadrupole interaction potential between a 
Ceo monomer at n and a surrounding alkali dumbbell at 
UA is given by 



UQ^{n,^{i7);nA,H'^A)) 



a,/3 



QaQl3 



A{n,a,'4!{v);nA,(3,Tl'{vA)) 



(4.30) 



where 



A(n, a, Tp{v);nA, /?, ipivA)) 

= \D{a,^{i^)) + X{n - ua) - dA{l3,^PiiyA))\ (4.31) 

is the distance between a charge Qa on the monomer and 
a charge Qp on the alkali dumbbell. The position of the 
charge /3 on the dumbbell ua in the crystal reads 



XinA,p)^X{nA) + dAiP,^ii^A)), 



(4.32) 



where dA is the position vector of the charge with respect 
to the center of the dumbbell. Explicitly we have 



dAi(3,i>{vA)) = ±dA {cos[V'(j7A)]e;- 



-sin[?/;(j74)]ec}, 

(4.33) 



where -I- and — refer to the two quadrupolar charges. 

The potential ( 4.3C ) is expanded in SARFs. After sum- 
mation over the polymer chains n an d the surrounding 
alkalis tIa we find (compare with Eq. (2.5a)): 



V^' ^ E E ^('^I'^!'^^l'^^)'^2(l7)s2( 



(4.34) 



where 



2ir 



2ir 



X{ni,i!]nAi,VA) ^ — I dip{iy) / d'ipipA) 
1- Jo Jo 

X U^'^in,4,{i;);nA,^{i^A))S2{4>{iy))s2iiyA)). 
Going over to Fourier space we have 



where 
A(gj 



AXh cos 



2 



2Xr COS 



(4.35) 



(4.36) 



(4.37) 



Here we have restricted ourselves to the six nearest neigh- 
bors Ha of a given monomer H. With the values dA of 
Table IV and with the same value qa = 0.12 for the three 



compounds, we have calculated the interaction energies 
Af, and Ac quoted in Table IV. The quantity |A(gj_)| is 
maximum at q±_ = qr in contradistinction with the di- 
rect RR interaction JQQ{qi_), Eq. ( 2.12] ). The intraionic 
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restoring forces of the electronic shells of the cations are 
described by a sum of single particle energy terms 



TABLE VI. Interactions at (fx = 9z, units K. 



'UQ±>2{q±), (4.38) 



with QA > 0. The self-energy gA is inversely propor- 
tional to the quadrupolar electronic polarizability and 
hence gcs < gm < gK (see Table IV). These con- 
cepts are inspired from the shell model of lattice dynam- 
ics where anisotropic electronic polarizabilities have been 
introduced The direct interchain coupling of alkali 
quadrupolcs is numerically small and will be neglected. 
We now consider the sum 



(4.39) 



where the subscript ep refers to the quadrupolar elec- 
tronic polarizability of the cations. We assume that the 
induced cation quadrupolcs follow adiabatically the mo- 
tion of the CgQ chains. For a given configuration {S2{q_L)} 
of the latter, we minimize with respect to S2{qi_) 

and find 



1 A(gl) 

2 9A 



S2{q±)- 



(4.40) 



Substitution into leads to the cation quadrupolar 

electronic polarizability mediated rotational interaction 



V. 



where 



RR 



N, Jep{q±)Sl{q±)S2{q±), (4.41) 



Jep{q±) = 77 



1 A2(gl) 



9A 



(4.42) 



This interaction is always attractive and maximum in 
absolute value at — qr (as is the case for the lattice 
mediated interaction V^^"). Both V^^ and V^J^ lead to 
a condensation of S2{q\_) at (fr and hence to the ferroro- 
tational structure /2/m: 



{S2{qA^) = 



(4.43) 



where a is the order parameter amplitude. 

In the following section we will discuss the competition 
between the direct rotational interaction of Cgo chains 
and the indirect, alkali mediated interaction. 



V. STABILITY OF STRUCTURES 

The orientational quadrupolar interaction be- 
tween Ceo chains is the sum of the direct interchain 
potential Vqq, Eq. ( 2.11 ) and of the indirect, alkali- 

" Eq. 



media ted potentials Vjy\ Eq. (4.27), and 



(4.41). Hence we write 



iRR 









— Jac{qz) - 


Jep{qz) 


J{qz) 


KCeo 


-195.84 


-184.70 





-49.26 


-429.80 


RbCeo 


-190.21 


-176.40 





-185.80 


-552.41 


CsCeo 


— 188.54 


— 167.88 





-227.64 


—584.05 


Ceo 


/ 


-121.86 


/ 


/ 


-121.86 




yRR 


_ yRR , 

- ^QQ + 


yRR 1 yRR 
* ac ~^ *^ ep 




(5.1) 



or equivalently 

rRR _ Nl 



Y,Jiq±)sl{q±)S2 q-L, 



with 



J{q±) = JQQ{q±) - Jaciq±) - Jepiq±)- 



(5.2) 



(5.3) 



Here the direct interaction JQQ{qi_), due to van der 
Waals and Coulomb quadrupole-quadrupole potentials, 
i s give n by En . ( 2.12| ), with Ja and Jt defined in Eqs. 
(|3.11a ), ( B.llbf ). The acoustic lattice dis placem ents me- 
diated interaction Jac{q±) is given by Eq. (4.28) while the 
electronic polarizability mediated interaction Jep{q±) is 
defined by Eq. ( 1.42| ). We have already shown that the 
direct interaction JQQ{q±) becomes maximum and at- 
tractive at the Brillouin zone boundary q± = qz, and 
hence favors the antiferrorotational structure Pmnn. On 
the other hand the alkali mediated interactions Jaciq±) 
and Jep{q±) both are attractive and maximum at the 
Brillouin zone center q± = qr, and hence favor the ferro- 
rotational structure /2/m. The strength of the indirect 
interactions depends on the specific nature of the alkali 
ions. 



We first consider J{qr)- From Eqs. ( |3.12D , ( [4.29[ ) and 
(4.42) we find the numerical values quoted in Table 0. 

Next we calculate J{qz)- The results are quoted in 
Table VI. We have taken into account that there is no 
coupling to acoustic phonons at qz- In Fig. 3 we have 
plotted the total rotational interaction J(q±) as a func- 
tion of along the line qr — qz- In Eq. ( p^ ) we have 
used for Jac{q\_) the interpolation expression 



Jac{q± = (0,9c)) 



8(z;,^J^sin^(gcf) 
q^Ci^abc 



(5.4) 



which coincides with Eq. (4.29) at qr and vanishes at 
qz- While for KCgo, with small polarizability of the K+ 
ion, the direct interchain interaction Jqgiqz) dominates 
and hence leads to Pmnn, for RbCgo and even more for 
CsCeo, the alkali mediated interaction Jep{qr) is most im- 
portant and leads to the structure /2/m. This structure 
is monoclinic and indeed one has measured small devi- 
ations of the (6, c) angle from 90 ° ||l3|,0 . The present 
theory accounts for such shears. We rewrite expression 
(|^ for (f_L = (0,gc) as 
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0,5 1.0 1.S 2.0 2.S 3.0 



q.c/2 

FIG. 3. Interchain energy J{q±), units K. 



We observe that 



lim iqc^Sr, (0,0, gc) = VniiVec,,, 
9c: — >o 



(5.5) 



(5.6) 



where e^^ are the homogeneous shears, and — NiNc is 
the nu mbe r of unit cells. Taking then the limit q,^ — > 
in Eq. (|^), we get by using Eq. ( [l.43| ) 



44 



(5.7) 



where a is the order parameter amplitude. We see that 
ferrorotational order induces 77 shears, in accordance 
with the monoclinic space group I2/m. With the numer- 
ical values of from Table the lattice parameters of 
Table ^ and assuming cr = 1, we obtain for e^,, the values 
of 0.476 ° and 0.483 ° for RbCeo arid CsCeo respectively. 
The deviations from 90 ° of the (5o,Co) angle measured 
experimentally are 0.316 ° ^ and 0.180 ° for Rb- 
and CsCgo respectively. 



VI. CONCLUSIONS 

In comparing the numerical values of the various 
contributions of the total interaction energies in Ta- 
bles and ^ we come to the conclusion that the 
quadrupolar electronic polarizability mediated interac- 
tion Jep((fj_) and the direct quadrupolar interaction 
jQQiq±) = ^QgiQi-) + Jqq{^1-) are the determining fac- 
tors for the structures Pmnn or /2/m, while the TR 
mediated interaction JaciQ±) plays quantitatively a neg- 
ligible role. The direct quadrupolar interaction, both of 
van der Waals and Coulomb type, is always attractive 
and maximum at q± = qz and hence favors the anti- 
ferrorotational structure, in agreement with the space 
group Pmnn for pressure polymerized Cgo p3|,p4| and 



for polymerized KCeo Notice that the values of 

Jqqiqz) for all three ACgo compounds are rather close 
together, and hence there is no way to explain the dif- 
ferent structure 12 /m for RbCeo ll|,|ll and CsCeo ^ 
by the direct interchain coupling via Cgo monomers, the 
only difference between the three compounds being the 
lattice constants. The role of the alkalis has to be more 
specific than just providing different lattice spacings. 

Since the values of JaciQr) are rather small and sim- 
ilar for the three compounds (Table 0), the only deci- 
sive interaction allowing to discriminate between KCgo 
and RbCgo, CsCeo, is the quadrupolar polarizability of 
the cations, introduced in Sect. IV. Here we do not in- 



voke the quadrupolar polarizability of the inner shells 
with small ionic radii but the quadrupole formation due 
to excited c?-states which carry a residual electronic 
charge. The average radii dA of these states, calculated 
in Appendix |^ and quoted in Table IV, are consider- 
ably larger than the ionic radii of K+ , Rb"*" , Cs"*" in con- 
ventional compounds like the alkali halides or like the 
ionic crystals with small molecular ions (alkali-cyanides, 
-nitrites). Within our view the large interstitial space 
available for the "octahedral" alkalis in ACgo compounds 
is a prerequisite for the existence and partial occupancy 
of these d-states in the solid. In that respect the alkali- 
fuUerides are unique ionic solids, and we expect that the 
excited rf-states are relevant for an understanding of the 
electromagnetic properties [H . Notice that in the present 
approach RbCgo and CsCgo are similar because their av- 
erage radii of excited Ad- and 5d-states respectively are 
close, while KCeo with an excited Sd-state has a consid- 
erably smaller value of . A same trend is also present 
in the atomic radii of the transition metals Sc, Y and the 
inner transition metal La, with 3c?^ , id^ and 5d^ electrons 
respectively, with atomic radii 1.64 A, 1.82 A and 1.88 A 
respectively. 

Within the present work we have assumed a model 
of rigid polymer chains which describes the structures 
Pmnn and /2/m. The assumption of rigid chains has to 
be abandoned if one wants to take into account a mod- 
ulation of electronic and structural properties along the 
orthorhombic a axis. In fact the interaction of electric 
quadrupoles between Ceo monomers along a same chain 
immediately suggests an antiferro-orientional structure 
along the chains. However such an extension is not yet 
sufficient to understand the recently discovered super- 
structure {a + c,b,a — c) below T = 50 K in KCeo • 
Here again the alkalis seem to play a specific role. 
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APPENDIX A: 

In the cubic phase the elastic part of the free energy 
per formula unit ACeo reads 

Uel = Y [cil {el^ + + el^) + 2ci2 {exx^yy + ^yy'^zz 



+ ^zz^xx) + 4c44 {e-ly + e^^ + ( 



(Al) 



Here Vc — a;;/4 is the volume of the primitive unit cell, Cc 
is the cubic lattice constant, are the elastic constants 
(we use the Voigt notation) and eij are the lattice defor- 
mations. Performing a rotation of the cubic system of 
axes {x,y,z) to an orthorhombic system (^,?7, C)i where 
^ corresponds to the cubic [110] direction, 77 to [TlO] and 
C to [001], we get 



Here the orthorhombic elastic constants c°,- are related to 

''J 

the cubic elastic constants by c\-^ = (cn + C12 + 2c44)/2, 

C12 = (cil-|-Ci2-2c44)/2, egg = Cii, C^3 = C12, cjg = C12, 

"-ee ~ ('-11 ~ ci2)/2, C55 = C44. Polymerization acts as a 
compression along [110], i. e. along ^. The corresponding 
surface forces are described by the boundary condition 

i5i 



(A3) 



where aik is the elastic stress tensor (indices i, k run over 
the orthorhombic axes ^,77,C); ^ is the normal to the 
surface. Using 



1 dUei 
Vc deik ' 



(A4) 



and observing that F = (— X, 0, 0), where if > is the 
strength of the compressional force, we obtain 



Since there are no lateral forces, = 
which leads to 

Cnew + c?2e« + c?3ecC = 0' 



(A5a) 
and <jQQ — 0, 

(A5b) 
(A5c) 



and a^n = 0, a^n = 0, tr^ij = lead to e^^ = 0, € (^ri = , 
e^^ — 0. Solving the system of equations ( A5a) -(|A5c| ) 
we ge t the deformations that are quoted in Eqs. ( p.laj )- 



APPENDIX B: 



TABLE VII. Born-Mayer-van der Waals potential con- 
stants for the interactions between the various types of ICs (a 
stands for "atom" , db for "double bond IC" and sb for "single 
bond IC"). The factors \ and arise from the threefold mul- 
tiplicity of the double bond ICs. For comparison with earlier 
work [34,50], Ci is expressed in units 3.745 x 10^ K and B in 
units 3.054 x 10^ K A^ C2 is given in A"\ 





a ~ a 


db ~ db sb — sb 


a — db 


a — sb 


db — sb 


Ci 


0.864 


i X 0.259 0.158 


i X 0.169 


0.169 





C2 


3.6 


3.2 3.6 


3.4 


3.6 





B 


2.1 















1. Born-Mayer and van der Waals Interactions 
between Ceo Monomers 

A Ceo monomer in ACeo (A=:K, Rb, Cs) is described 
by a rigid cluster of interaction centers (ICs). The ICs 
refer to pair potentials between monomers on separate 
chains. In addition to the 60 carbon atoms, three ICs 
per double bond (at positions L/2 and ±L/4, where L 
is the bond length) and one IC per single bond (located 
at the bond center) are considered. Following we 
consider van der Waals potentials of type ( p.q ), where 
the potential parameters Ci, C2 and B depend on the 
type of ICs. Numerical values are given in Table VII. 
The comparison of experimental X-ray diffuse scattering 
results with theoretical predictions based on various in- 
termolecular potentials, shows that such a model has its 
merits ] ^^ . However the transition temperature for 
the phase change FniSm — > Pa3 obtained by poten- 
tials used in [|^ is much too large |42 4^ in comparison 
with the experimental value Tc ~ 250 K. It is then nec- 
essary to adjust the potential parameters. A potential 



is considered satisfactory if it reproduces the transition 
temperature and the experimental crystal field [Q in 
the orientationally disordered phase of Ceo- The poten- 
tial parameters of Table VII give Tc = 252 K and crystal 
field expansion parameters wq = 346.3, wio = —84.1, 
u'12,1 = —82.9, wi2,2 = —378.6 for the disordered phase 
of sohd Cqq. 



2. Born-Mayer and van der Waals Interactions 
between a Ceo Monomer and an Alkali Atom 

Here we let the alkali atom interact with the 60 car- 
bon atoms of the Ceo monom er. We take again the pair 
potential of the form of Eq. ( ^.5[ ) . Numerical values for 
the potential constants are quoted in Table VIII. These 
values are derived from the homoatomic interaction po- 
tential constants: for C— C the a — a values of Table VH 
were used, for K— K and Rb— Rb the values of Ref. |45| ] 
were taken, and for Cs— Cs an extrapolation of the K— K 
and Rb— Rb values based on the ratios of the ionic radii 
of K+, Rb+ and Cs"*" was made. 
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TABLE VIII. Born-Mayer-van der Waals potential con- 
stants for the interaction between a carbon atom and an alkali 



atom. 




units 


C-K 


C-Rb 


C-Cs 


Ci 


10' K 


2.43 


3.59 


4.96 


Ca 


A-i 


3.28 


3.29 


3.32 


B 


10^ K 


3.36 


5.26 


8.23 



APPENDIX C: 

To calculate the portion of the charge located inside 
the interstitial sphere precisely, one should perform an 
electron band structure calculation of ACgo- However, 
some conclusions can be made on the basis of a general 
consideration. In the tight-binding approximation we ex- 
pect that the electron wave function of a band electron 
at an alkali site n, ipAik, a)\fi, is expanded in terms of 
local s and d functions, i.e. 



ipA{k,a) 



1 



Jk-X{n) 



N 

2 

E 



7oo(fc, a)Rs{r)Yo° 



l2rr.{k,a)Rd{r)Yr{f) 



(CI) 



where Rs{r) and Rd{r) are the radial parts of s and d 
states of the alkali, respectively. {N is the total number 
of sites in the crystal.) Here we are speaking about 4s 
and 3d electron states for K, 5s and Ad states for Rb, 
and 6s and Sd states of Cs. The s and d coefficients 
7i,„(fc,Q!) (i.e. 7oo(fc,a) and 72m(fc,a)) are found by 
solving a secular equation for the band electron with the 
wave vector k and the band number a. The inclusion 
of the d-states in ■4!A(k,a)\fi is logical since for a neu- 
tral K, Rb and Cs atom the d-shell corresponds to the 
first excited electron level. The resulting charge distri- 
bution inside the sphere at the site n can be computed 
by summing up the density of all extended states fc, a 
below the Fermi energy, E{k,a) < Ep. The charge at 
the alkali site is small but not zero and we believe that 
|7g(fc,Q;)P ~ 0.1, |7^(fc,Q;)p - 0.1. It is important that 
the wave function ipA{k, a) has nonzero matrix elements 
of quadrupolar charge density. To demonstrate this, we 
consider quadrupole components of the electronic den- 
sity associated with SARFs 5'A(r), where A refers either 
to the two components (fc = 1, 2) of Eg symmetry or to 
the three components (fc — 1,2,3) of symmetry, i.e. 
A = (Z = 2,r, fc), where T = Eg or Tzg. In Ref. ||6| it 
has been shown that the operator of quadrupolar density 
of conduction electrons reads 



k 



(C2) 

where li, I2 = or 2 (s and d states), and at , a^^ are 
creation and annihilation operators for one electron in 



the state {k,a). Here 

mi ,7712 

xc\{limi, l2m2). (C3) 
with the quadrupolar matrix elements 



CAihmi 



12^2) = j 



\n)SA{^)Yi"'^{n) dn. 



(C4) 



Some coefficients CA(^imi, 121^2) are different from zero 
for certain d — d and d — s transitions as has been 
shown in Ref. p6| within a more general context (see 
there Eq. (3.21) for s — d transitions, and Table || for 
d — d transitions with A = (T2g,k)). From the latter 
observation we conclude that if there is an appreciable 



weight of d-states in ■(/'^(fc, a)|^;, Eq. (CI), then the al- 
kalis acquire a quadrupole moment. (In order to obtain 
the operator of the quadrupolar moment for the ground 
state of A Cfio we put q = k — p = and a = P in 
Eqs. (C2) and (|C^).) This consideration gives a micro- 
scopic support for the concept of quadrupolar polariz- 
ability of alkalis introduced in Ref. |Q. It is clear that 
the quadrupolar moment of alkalis depends on the co- 
efficients 700(^,0;), l2m{k,a). In principle, the secular 
equation for the coefficients should also take into account 
the bilinear quadrupole-quadrupole interaction with the 
polymer chains of Cgo molecules. In practice, however, 
such calculation would be very difficult to implement. In- 
stead, following Ref. one can introduce phenomeno- 
logically quadrupolar moments of alkalis and then op- 
timize their interactions with the neighboring polymer 
chains. In this work we use a simple model with two 
point charges qA at distances idyi from the alkali center. 
Here dA is an average d-shell radius which is found as 



{r d 



r'^dr' . 



(C5) 



The radial functions TZd and the average radii {dx, dub, 
dcs) were calculated numerically. We have employed a 
relativistic atomic program for self-consistent-field cal- 
culations of K+, Rb"*" and Cs"*" with the local density 
approximation (LDA) of exchange according to Barth- 
Hedin Since the ionicity of alkalis is close to +1, 

TZd{r) corresponds to the virtual (almost empty) 3d3/2 
shell for K, to the virtual 4d3/2 level for Rb and to the 
virtual ^dz/2 level for Cs. The role of excited d-functions 
has been discussed by Murrel and by Niebel and Ven- 
ables ]49| ] in connection with the problem of explaining 
the observed crystal structure of the rare gas solids. 

On the other hand, as it has been shown in Ref. ^] 
there is electronic on-site interaction which is propor- 
tional to the square of the quadrupolar moment (see Eq. 
(3.14a) of Ref. |46|). This interaction gives rise to the in- 
traionic self-energy qa introduced in Eq. (6) of Ref. pO] 
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and Eq. (4.38) of the present work. Upon uniform ex- 
pansion or contraction of the electron density around an 
alkah, the strength of it scales as l/dA- 
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